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OPTIMIZING VENTILATORY ASSIST IN VENTILATOR-DEPENDENT 
PATIENTS USING ELECTRICAL ACTIVITY OF RESPIRATORY MUSCLES 



FIELD OF THE INVENTION 



The present invention relates to a method for determining an optimal 
level of ventilatory assist to a ventilator-dependent patient, using myoelectrical 
activity. 



10 BACKGROUND OF THE INVENTION 

Of the many factors that lead to respiratory muscle fatigue, the tension 
developed by the muscles [34] and the duration of a contraction [2] are 
particularly important. These two factors are contained in such indices as the 
tension-time index [3] and the pressure-time product [10, 20, 32,. 35]. 

15 Bellemare and Grassino [3] showed a direct inverse relationship between the 
endurance time of a fatiguing diaphragm contraction and the rate of decay of 
the ratio of the high to low spectral components (H/L) of the diaphragm 
electrical activity (EAdi), indicating that the two are expressing progressive 
failure to sustain load. The force exerted by the muscle has been shown to be 

20 directly related to the rate of decay of the power spectrum center frequency (or 
H/L ratio), and the level at which it plateaus [16, 21, 28]. Such shifts in the 
power spectrum reflect a slowing in the muscle action potential conduction 
velocity [28, 38, 39], and constitute an early indication that, at the cellular level, 
such breathing patterns cannot be maintained indefinitely [3]. 

25 

• Hyperinflation, which impairs the length-tension relationship of the 
respiratory muscles, i.e. the transformation of the neural activation into a 
mechanical output or pressure, reduces the capacity of the respiratory 
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muscles to generate pressure (neuromechanical uncoupling), unless EAdi is 
increased. Studies have shown that when inspiratory pressure, flow and duty 
cycle remain constant, increases in end-expiratory lung volume (EELV) 
promote reductions in endurance time [33, 44] and sustainable pressures [1 1]. 
5 In an animal model, Tzelepis et al [44] proposed that under such conditions, 
diaphragm shortening would require greater excitation to generate a given 
sub-maximum tension, and that this increased activation might account for the 
greater contractile muscle fatigability observed at shorter muscle lengths. 

10 The "optimal level" of partial ventilatory assist, with the aim to ensure 

adequate pulmonary ventilation while preserving inspiratory muscle function, is 
generally set on an empirical basis in the clinical setting. It has also been 
proposed that the optimum level could be determined from the lowest stable 
breathing frequency (f B ) achieved, i.e. without bradypnea or apnea. In the 

15 patients, this corresponded to 16.4 bpm (breaths per minute) and was 
associated with a tidal volume (V T ) of 1 1 -8 ml/kg. However, mechanical lung 
modeling in that study demonstrated that such a level of support actually 
resulted in a total unloading of the respiratory muscles. Others have defined 
the optimum level of partial ventilatory assist as that which produces the 

20 lowest transdiaphragmatic pressure (P d i) swings and found that this was 
associated with a f B of 19.7 bpm and a V T of 11.7 ml/kg. The pressure P dl in 
the study was used as a marker of inspiratory effort. Jubran et al [20] defined 
an upper bound inspiratory pressure-time product <125 cm H 2 0 s/min as the 
desirable level of inspiratory effort to be achieved during partial ventilatory 

25 assist. Although arbitrarily determined, the appropriateness of such a threshold 
was justified by the fact that it corresponded to a tension-time index (TT d i) that 
was well below that considered to signify impeding inspiratory muscle fatigue. 
The study found a high variability in pressure-time products between patients 
and demonstrated that a f B < 30 bpm and V T of 0.6 L were better determinants 

30 of an. optimal level of inspiratory effort during partial ventilatory assist. Based 
on these breathing pattern findings, it is likely that the level of respiratory 
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• muscle unloading provided by this method of optimizing partial ventilatory 
assist was lower than that of the above reviewed studies. Brochard et al [8] 
defined the optimal partial ventilatory assist level as the lowest level of 
ventilatory assist, which when implemented, maintained the highest level of 
5 diaphragmatic electrical activation without the occurrence of fatigue as 
evaluated via diaphragm EMG power spectrum analysis. Interestingly, such 
levels of partial ventilatory assist were associated with a f B of 20-27 bpm and 
V T of 8.0 ml/kg, values similar to those later reported by Jubran et al [20]. . 

10 

SUMMARY OF THE INVENTION 

In accordance with the present invention, there is provided a method 
for determining an optimal level of ventilatory assist to a ventilator-dependent 
15 patient, using myoelectrical activity. For that purpose and according to non- 
restrictive illustrative examples, the method uses one or a combination of the 
following} relations: 

- A first relation estimating the patient's diaphragmatic muscle force from the 
20 patient's EAdi; 

- A second relation estimating the patient's diaphragmatic muscle force from 
spectral changes during diaphragm contraction and involving a critical 
force level above which muscle fatigue starts to develop; 

25 

- A third relation, deduced from the first and second relations and relating 
the patient's diaphragmatic muscle force to (a) the critical force level above 
which muscle fatigue starts to develop, and (b) values related either to the 
patient's EAdi or the spectral changes during diaphragm contraction; and 

30 

- A fourth relation relating the transdiaphragmatic pressure to either (a) the 



\ 
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patient's EAdi or (b) the spectral changes during diaphragm contraction. 

The above relations will enable evaluation of many patient's respiratory 
features from measurements of values such as EAdi, gastric pressure, 
esophageal pressure, respiratory flow, etc. The evaluation of these patient's 
respiratory features will allow the medical professionals to better optimise the 
characteristics of assistance of a ventilator-dependent patient. 

The foregoing and other objects, advantages and features of the 
present invention will become more apparent upon reading of the following 
non-restrictive description of illustrative embodiments thereof, given by way of 
example only with reference to the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 
. In the appended drawings: 

Figure 1 is a schematic representation of a non limitative example of 
experimental setup. Left A catheter-mounted multiple array esophageal 
electrode used to measure diaphragm electrical activity (EAdi) with balloons 
mounted upon it for measurement of esophageal pressure (P es ) and gastric 
pressure (P ga ) was passed transnasally and positioned at the gastro- 
esophageal junction. Airflow was measured with a pneumotachograph and 
tidal volume (V T ) was obtained by integrating inspiratory flow. Right On-line 
display of the target transdiaphragmatic pressure P dj and the root-mean- 
square (RMS) of the EAdi. 

" Figure 2 illustrates examples of tracings of tidal volume (V T )3aphragm 
electrical activity (EAdi), transdiaphragmatic pressure (P di ), esophageal 
pressure (P es ), and gastric pressure (P ga ) obtained in one subject during 
"volume" maneuvers and "expulsive" maneuvers performed during the 
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hereinafter reported study. The "volume" maneuver consisted of an end- 
inspiratory hold at an increased lung volume, which resulted in the generation 
of a low P dl (left tracing), whereas the two expulsive maneuvers were 
performed at end-expiratory lung volume targeting a low P d j (middle tracing) 
5 and high P dj (right tracing). 

Figure 3 shows examples of bar graphs displaying drops in diaphragm 
power spectrum center frequency (CF dj ), targeted levels of P dj , pressure-time 
product (PTP dj ) and the associated electrical activation (EAdi) observed during 
10 the three maneuvers performed in the study. Bars are average values 
obtained for five subjects (± SD). 

: Figure 4 represents examples of graphs from one representative 
subject showing the center frequency (CF dI ), root-means-square (RMS) of the 
15 EAdi and the diaphragm pressure-time product (PTP d j) plotted over time during 
the volume maneuver (circles) and the two expulsive maneuvers at EELV, one 
targeting a low P di (squares) and the other a higher P di (triangles). 

20 DETAILED DESCRIPTION OF THE ILLUSTRATIVE EMBODIMENTS 

the study was conducted to determine in humans whether an 
increased EAdi, with neuromechanical uncoupling, promotes greater 
reductions in the center frequency of the EAdi signal (CF d i), when the 
diaphragm pressure-time product (PTP d i) is kept constant. An additional aim* of 
25 the study was to establish the extent to which PTP d | would need to be 
increased, in the presence of normal neuromechanical coupling, in order to 
reproduce the drop in CF d i observed with uncoupling. 

More specifically, the study evaluated whether increased diaphragm 
30 activation induced by an increased lung volume promotes increased drops in 
the center frequency (CF dl ) of the diaphragm's electrical activity (EAdi) when 
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pressure-time product (PTP dj ) is kept constant. Five healthy subjects 
performed runs of intermittent quasi-static diaphragmatic contractions with a 
fixed breathing pattern. In separate runs, subjects targeted transdiaphragmatic 
pressures (P dl ) by performing end-inspiratory holds at total lung capacity with 
5 the glottis open (neuromechanical uncoupling), and at end-expiratory lung 
volume by performing expulsive maneuvers (no neuromechanical uncoupling). 
Diaphragm activation and pressures were measured with an electrode array 
and with balloons mounted on an esophageal catheter, respectively. 
Reproducing a P di of ~31 cm H 2 0 during neuromechanical uncoupling (lung 
10 volume increased to 77.5% of the inspiratory capacity) increased EAdi from 
25% to 61% of maximum (P<0.001) and resulted in a 17% greater drop in 
center frequency (Fb=0.012). In order to reproduce, in the absence of 
neuromechanical uncoupling, the decrease in center frequency observed 
during neuromechanical uncoupling, a two- fold increase in P dl and PTP dj was 
15 required. It was concluded that a constant pressure-time product does not 
guarantee that the center frequency of the diaphragm remains stable when 
activation is increased. 

METHODS 

20 

Subjects 

Five healthy volunteers (1 female, 4 males) with a mean age of 40.6 ± 
8.0 years participated in the study. The study was approved by the Scientific 
25 and Ethical Committees of Sainte-Justine's Hospital and all subjects gave their 
informed consent. 

Experimental protocol 

30 While seated in an upright chair, and facing a computer monitor 1 

(Figure 1), each subject 2 performed repeated maximal inspirations to total 
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lung capacity (TLC) in order to obtain three reproducible voluntary maximum 
EAdi values. Subjects 2 were subsequently asked to perform intermittent, 
near-isometric diaphragmatic contractions of 10 seconds duration, separated 
by 5 seconds relaxation periods in which free breathing was allowed. With 
5 visual feedback of transdiaphragmatic pressure (P d} ), during two runs a low 
level of P di was targeted, while a higher level was targeted during the third run. 
The duty cycle was imposed by a sound signal, and each run lasted until a 
plateau in CF dl was reached, or until the subject was no longer able to 
maintain the target P di . 

10 

In order to obtain two different EAdi levels for the same target P di , 
subjects were instructed to perform two different maneuvers: 

1 . Volume maneuver: Subjects inspired close to their TLC and produced a 
15 given level of P dj (Figure 2; left tracing). The P di was maintained at this lung 

volume with the glottis open. 

2. Expulsive maneuver: Subjects performed expulsive maneuvers in order to 
generate a target P di . All expulsive maneuvers were performed at end- 

20 expiratory lung volume (EELV) (Figure 2; middle and right tracings). 

After initially performing a volume maneuver run, subjects then 
performed two expulsive maneuver periods. One expulsive maneuver run 
targeted a P dj similar to that observed during the volume maneuver but 
25 requiring less EAdi, and the other targeted an increased P dl (to reproduce the 
CF dl observed during the volume maneuver). The volume maneuver was 
subsequently repeated once for retest purposes. Subjects rested for 20 
minutes between all runs. 

30 Instrumentation 
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Using the set-up of Figure 1 , airflow and tidal volume were measured 
by the computer 4 through a pneumotachograph 3, EAdi was measured by the 
computer 4 through a linear array of electrodes such as 5 mounted on an 
esophageal catheter 6 inserted through a patient's nostril (or patient's mouth) 
5 until the electrodes 5 are positioned in the gastro-esophageal junction 10 of 
the patient's diaphragm 7, esophageal (P es ) and gastric (P ga ) pressures were 
measured by the computer 4 through gastric 8 and esophageal 9 balloons 
mounted on the cathether 6 on opposite sides of the array of electrodes 5. 
Transdiaphragmatic pressure (P dl ) was obtained by the computer 4 by 
10 subtracting P es from P ga . 

On-line automatic processing of EAdi 

EAdi (calculated as the root-mean-square (RMS)) was acquired, 
15 processed and displayed on-line using a standardized methodology [4, 36, 41], 
CFdi was evaluated for signal quality using established indices and criteria in 
accordance with Sinderby et al [40]. To avoid influence of power spectral shifts 
on EAdi signal strength, RMS was calculated on the spectral moment of order 
1 (Ml) which is insensitive to conduction velocity [6]. For more extensive 
20 review see Aldrich et al [1 ]. 

Off-line signal analysis 

Inspiratory duration (T,) f total breath duration (T, ot ), and breathing 
25 frequency (f B ), EAdi and pressures were determined using the P di signal. The 
PTP d i was obtained by multiplying the area subtended by the P di signal by the 
respiratory frequency. EAdi signal amplitude was expressed as a percentage 
of the voluntary maximum EAdi (obtained from TLC maneuvers) [37]. 
Variables were compared between each of the maneuvers performed using 
30 one-way repeated measurements analysis of varianc" (ANOVA) and post hoc 
contrasts of significant effects were performed using the Student-Newman- 
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Keuls test. Test-retest reliability of the P di , EAdi and CF d i values obtained 
during the two volume maneuvers performed was evaluated by calculating the 
interclass correlation coefficient (ICC). 

5 RESULTS 

The subjects were able to perform all maneuvers and maintain the 
imposed duty cycle (P=0.93; one -way ANOVA) during all protocols (Table 1). 

10 As shown in Table 1 and Figure 3, subjects were able to achieve and 

maintain similar target P di levels during the volume maneuver (high lung 
volume)' and the low-pressure expulsive maneuver (EELV) (P=0.99). During 
the volume maneuver, subjects inspired to an average of 77.5±11.1% of their 
inspiratory capacity. In order to generate the same PTP dj at different lung 
15 volumes, the volume maneuver (neuromechanical uncoupling) required an 
EAdi of 60±8% of maximum compared to 25±8% for the expulsive low- 
pressure maneuver (at EELV). As shown in Table 2 and Figure 3, despite a 
matching of the PTP di , the volume maneuver promoted a 17% larger drop in 
the CF di than the expulsive low-pressure maneuver (EELV). Figure 4 shows a 
20 time tracing from one representative subject, where CF dj declines more rapidly 
and to a greater extent during the volume maneuver (circles), which required a 
high EAdi for a similar PTP dl compared to the expulsive low-pressure 
maneuver. 

In order to produce a similar drop in CF di during the expulsive 
maneuver (EELV) as was observed during the volume maneuver, more than a 
two-fold increase in the target P dl (P=0.015) was required. This was associated 
with an increase in EAdi from 25±8 % to 44±9% of maximum J^O.009) 
(Figure 3). As can be seen in Figure 4, the rate of decline of the CF di was 
similar for the volume maneuver (circles) and the expulsive high-pressure 
maneuvers (triangles). 



25 



30 



x 
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Presented in Table 3 are the values of P di , EAdi and CF di for the test- 
retest of the volume maneuver. During the retest, subjects successfully 
targeted a P di that was similar to that generated during the initial volume 
5 maneuver (ICC=0.95). EAdi was also similar (ICC=0.93) as was the drop in 
CFdi (ICC=0.98). 

DISCUSSION 

10 

. The study evaluated intermittent static contractions maintained at two 
different lung volumes, in order to examine the effect of altered 
neuromechanical coupling and increased diaphragm electrical activation, on 
diaphragm sarcolemma excitability, assessed by changes in CF di . It was found 
15 that, for a given targeted PTP di , the drop in CF di was greater when EAdi was 
increased by neuromechanical uncoupling, suggesting that the level of muscle 
activation influences CF di . 



Studies on the canine diaphragm have demonstrated that changes in 
20 CF di are associated with changes in the mean action potential conduction 
velocity (APCV) [38], confirming previous mathematical models [29]. During 
muscle contractions, both CFdi and muscle fiber APCV depend to a smaller 
extent on the cable properties of the fiber, [38, 39], and to a larger extent on 
the muscle membrane excitability [17, 18, 29, 39]. The excitability of the 
25 muscle fiber membrane is dependent on the transmembrane gradient of 
potassium, and with increased muscle activation, efflux of potassium 
increases. In order to defend the extracellular potassium concentration (and 
hence, the membrane potential), the cell depends on the re-uptake of 
potassium, e.g. via the ATP (Adenosine Triphosphate) dependent 
30 sodium/potassium pump [12], and washout via the blood circulation [25], i.e. 
diffusion of potassium from the extra-cellular space into- the blood stream. 
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Regardless if blood flow is reduced [23, 31, 42], or the muscle activation is 
increased, as in the present work, the muscles electrical activity will indicate 
reduced, membrane excitability, by shifts in the power spectrum toward lower 
frequencies. CF di can also be affected by factors such as motor unit territory, 
5 number of fibers in the motor unit, dispersion in arrival times of the single 
contributions in the motor unit signal, dispersion in action potential conduction 
velocities between motor units that can cause the EAdi power spectrum to shift 
[4, 29]. However, given that these influences are minor in healthy muscles [30] 
and given that the test situation did not allow for much variability in the 
10 contractile pattern, it is unlikely that these influences had more than a minor 
impact on the results. 

In the study, a constant P dl - was targeted with a constant duty cycle at 
two different lung volumes, and it was therefore assumed that P di hindrance to 
15 blood flow under those conditions remained relatively similar at the different 
muscle lengths [19]. However, in order to achieve the same target P di at an 
increased lung volume, EAdi was increased, which represents an increase in 
energy demand/consumption as well as increased metabolic output (e.g. 
potassium efflux) from the cell. As can be seen in Figure 3, the rate of decline 

20 of CF di at increased lung volume was significantly higher than that observed 
when the same pressure was targeted at FRC with lower EAdi. Vitro studies 
have also demonstrated that increased activation (i.e. demand), accomplished 
by increasing stimulation frequency of a muscle shortened to 70% of its 
optimum length, in order to obtain the same tension generated at optimum 

25 length, resulted in an increased fatigue in the shortened muscle [14]. The 
current study therefore demonstrates that the higher diaphragm activation 
required for generating the target P di at an increased lung volume 
(neuromechanical uncoupling) influences the rate/extent to which CF di decays. 
Further theoretical evidence for the impact of neuromechanical uncoupling on 

30 the CF di is provided in the Appendix. 
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In the absence of neuromechanical uncoupling, an increase in P d j is 
always associated with an increase in EAdi. In the study, doubling of P dj at the 
same lung volume (i.e. FRC) was associated with an increase in EAdi from 
25% to 44% of maximum. Beck et al [6] showed that EAdi (in absolute values) 
5 is closely related to P dj , such that activation increases (i.e. energy demand 
increases) when pressure increases (i.e. energy supply decreases). However, 
this relationship is altered when muscle length is changed. In such a 
circumstance, P di continues to reflect EAdi only when the P di is normalized to 
the maximum P dl obtained at each corresponding lung volumes [6]. It was 

10 previously shown that when the same EAdi is targeted at different lung 
volumes, the higher resulting P di generated at FRC promotes a greater drop in 
CF dj than does the lower pressure produced at the higher lung volume [42]. 
Such results indicate that for a given neural activation, an increase in force or 
P di reduces diaphragm excitability. Consequently, the use of the TT dl and 

1 5 pressure-time product as indices for predicting changes in the excitability of 
the diaphragm sarcolemma (as reflected by CF dl ) is limited to conditions of 
constant neuromechanical coupling, where the diaphragm force generating 
capacity remains unaltered. 

20 Consistent with previous studies [3, 16, 21, 28], doubling of the target 

P di at FRC in the present study increased the rate of decline of the CF di as well 
as the level to which it declined (Figures 3 and 4). This is partially explained by 
the increase in EAdi, as discussed above. However, it is also partially 
explained by the fact that: 
25 i) diaphragm contractions with a high P dl tend to hinder blood flow (i.e. 

energy supply) relatively more than contractions producing a low P dl 
[19], and 

ii) impaired blood flow to a muscle has the propensity to promote shifts 
in the electromyographic power spectrum toward lower frequencies 
30 [22,30]. 
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Methodological and technical aspects 

In the study the contraction and relaxation periods were maintained at a 
fixed duration and therefore any potential influence of duty cycle on muscle 
5 function [2, 22] was controlled for. It must be emphasized that accurate 
physiological measurement of CF d i depends on being able to control for: 

(a) changes in muscle-to-electrode distance, 

(b) electrode positioning with respect to the muscle fiber direction and 
location, 

10 (c) electrode configuration, 

(d) signal to noise ratio, 

(e) influence of cross-talk from other muscles (including the heart and the 
esophagus), 

(f) electrode movement-induced artifacts [7, 36, 38, 39, 40]. 

15 

In the study, the technology used to measure EAdi power spectrum 
includes all means and devices necessary to minimize these influences [1 , 36, 
40]. The findings that evoked muscle action potentials are influenced by 
changes in lung volume [5, 15] have contributed to the assumption, of a 

20 potential inherent inaccuracy of measured EAdi amplitude [5, 15] and CF di [5]. 
However, during mild voluntary muscle contractions that do not alter 
diaphragm membrane excitability, it has been shown that chest wall 
configuration/lung volume and changes in muscle length do not have an effect 
on EAdi and CF d j [5, 6, 7, 17, 39]. Therefore the above-discussed effect of 

25 chest wall configuration/lung volume likely did not have an impact qn the 
results. 

Another factor that could have influenced our results was the difference 
in the partitioning of esophageal and gastric pressures for the same P d i during 
30 the different maneuvers. In a previous study [42], where subjects targeted the 
same EAdi at high and low lung volumes, greater decreases in CF dj were 
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consistently observed at EELV (higher Pdi), regardless of whether subjects 
performed an expulsive (i.e. Pdi generated mainly by gastric pressure) or a 
Mueller maneuver (i.e. Pdi generated mainly by esophageal pressure) at EELV 
[42], In a pilot trial to that study (unpublished observations), it was found that 
5 diaphragm contractions generating identical Pdi, duty cycle and EAdi, 
produced the same trajectory of decrease in CFdi, whether subjects performed 
expulsive or Mueller maneuvers. Therefore, it is not believed that differences 
in the partitioning of the esophageal and gastric pressures during the volume 
and pressure maneuvers in the current study had an effect on the outcomes 
10 observed. 

Clinical implications 

The results of the present study have direct implications to patients 
15 being weaned from mechanical ventilation. It is well known that patients 
undergoing a weaning trial may demonstrate dynamic changes in EELV 
(dynamic hyperinflation) [43], which similar to the study would alter the 
neuromechanical coupling of the diaphragm. In order to compensate for this 
uncoupling (i.e. maintain the same P di ), the patient would need to increase 
20 diaphragm activation. The combination of an increased activation, with an 
elevated P di would, according to the present study, lead to decreased 
diaphragm CF di (excitability), and possibly an increased respiratory effort 
sensation [42]. Shifts in the H/L ratio of the power spectrum of the diaphragm 
electrical activity have been reported in patients with respiratory failure in 
25 whom ventilatory assistance is removed [8, 13]. However, given that 
diaphragm weakness is prevalent in mechanically ventilated patients [24], and 
it remains to be determined at what combined levels of EAdi and P dj would 
affect.CF d i. 

30 Conclusion 
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The study shows that diaphragm activation is an important determinant 
of diaphragm membrane excitability, and changes in CF di . Furthermore it 
shows that the pressure-time product and tension-time index cannot be 
considered as valid reflections of diaphragm energy consumption and/or 
sarcolemma excitability when neuromechanical coupling is altered. 

Appendix 

With data from the above investigation, the diaphragmatic muscle force 
can be estimated from myoelectric measurements in two ways. The first way is 
expressed with the equation: 

F = ixEAdi (1) 

where Fis the diaphragmatic muscle force, fi is a proportionality constant, and 
EAdi is. a measure of the myoelectric signal strength of the diaphragm.. Here 
the square root of the first power spectral moment is used since it represents 
the signal strength, which has been compensated for the influence of changes 
in the propagation velocity of the myoelectric action potentials [29], 

The second way to estimate the diaphragmatic muscle force is from the 
spectral changes during diaphragm contraction. For a forceful periodic muscle 
loading, the initial center frequency CF 0 decreases to a final plateau value CF TO 
according to the equation [26]: 

CF„ - CF 0 (X-K)T D /[(1-K)T D +KT R ] (2) 

Here, /r is the duty cycle, Le. the inspiration time in relation to the'total time 
period, and T R is the CF d i recovery time constant pertaining . to an 
approximately exponential time curve which is rather independent of the 
muscle force [9]. The symbol T D denotes the time constant for the decrease in 



\ 
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CFdi, which is related to the muscle force as [27]: 

T D =t)KF-F c ) (3) 

5 In this equation rj is a proportionality constant and F c is a critical force level 
above which muscle fatigue starts to develop. Equation (2) is rearranged to 
obtain the experimentally determinable quantity: 

Q = t r IT D = [d - k)Ik ] [(CF 0 - CF„)/CF„] (4) 

Equations (3) and (4) then give: 

F = F C +Qr]/T R (5) 

15 Making equal the two force estimates of equations (1) and (5) the following 
relation is found: 

ocEAdi -p-Q = o (6) 

20 where 

a = i£T R lr\ (7) 

and 

25 

P=F c T K /rj (8) 

Relation (6) represents a set of three equations (the three experimental 
conditions) with two unknowns. A fitting procedure with data from Table 4 (with 
30 simultaneous minimization of the relative errors in EAdi and Q) gives the 
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values a = 0.00417 and ft = 0.0419 with a relative fitting error of 0.24. 



With a and (3 known, the experimental values of the force can be 
expressed in relation to the critical level for onset of deterioration of cell 
5 excitability. The two ways to describe this are obtained by rearrangement of 
equations (1) and (7), and equation (5), respectively, which gives: 

0 7 = (F/F C ), = aEAdi/J3 (9) 

10 and 

*,-<F/F c ) ff + (10 ) 

These quantities have been determined and are listed in Table 4 together with 
1 5 their mean values </> m . 

The observed transdiaphragmatic pressure (P d j) is assumed to be 
related to the muscle force as: 

20 P dl = FG (11) 

where G is a geometrical factor taking into account that the diaphragm muscle 
changes its shape with the inspired volume. This factor is thus assumed to be 
the same during the expulsive maneuvers with low or high P dl generation 
25 performed at end-expiratory lung volume. As with the force relations, the 
transdiaphragmatic pressure can be expressed in two ways, relating to the 
myoelectric signal strength and to the fatigue induced spectral changes. 
Combining equations (1), (5), and (1 1) we obtain: 

30 Pdi = fi EAdi G (12) 
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and 

Pdi = {F c +Qr}IT R )G (13) 

5 which is further developed (with equations (7) and (8)) into the two 
expressions: 

Yr^iGri/T^^Pdi/ia EAdi) (14) 

and 

10 

Y tt =iGri/T R ) ff =Pdi/(P + Q) 

(15) 

Numerical values, calculated for the two expressions, are given in Table 4 
together with their mean values y m . 

15 

From the results, listed in Table 4, it can be concluded that the forces in 
relation to the critical level are approximately the same during the volume 
maneuver and the high pressure expulsive maneuver, which is also reflected 
in their deterioration of cell excitability, expressed by the factor Q. During all 

20 conditions the forces are above the critical level as seen in the factor <p m . The 
geometrical dependence, expressed in the factor y m , is obviously the same 
during low pressure expulsive maneuver and high pressure expulsive 
maneuver, but is much less during the volume maneuver. The ratio between 
the y values in the volume maneuver and the expulsive maneuver is 0.41 . 

25 Since the rj values and the T R values are expected to be independent of the 
maneuvers, this means that also the G factors have the same ratio. This 
indicates a much lower efficiency to convert force into pressure during the 
volume maneuver. The tension time index, taking into account the timing and 
the pressure, is thus not sufficient to describe the complexity of fatigue 
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development. At least it has to be modified with a volume dependent 
correction factor. Better, though, are methods reflecting the deterioration of cell 
excitability and not the mechanical result of the contraction. 



5 Electromyographic and mechanical methods to detect muscle fatigue 

Based on the above results, methods to determine critical levels of 
muscle failure during periodic loading (such as respiration) are described. A 
number of equations relating certain physiological variables to each other are 
1 0 needed and they will be derived prior to the description of the methods. 

Periodic muscle load characteristics 

Consider a periodic muscle loading, such as the respiratory work, in 
15 which repeated muscle contractions alternate with muscle relaxations. The 
situation is characterized by the time period T 0 and its parts: the duration of 
each contraction Ti and the duration of relaxation T 2 where: 

To = Ti + T 2 (16) 

20 

In order to simplify the equations, the duty cycle k is introduced as: 
: k = Tt /T 0 (17) 
25 The mean force developed during the time interval Ti is denoted F. 

Myoelectric changes due to fatigue 

Isometric fatiguing contractions cause the myoelectric center frequency 
30 CF to decrease exponentially from its resting value CF 0 with a time constant 
T F . During recovery the center frequency returns gradually to its normal value 
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following an approximalety exponential course, described by the recovery time 
constant T R . It should be observed that many other characteristics of the 
power spectrum of the myoelectric signal exhibit the same dependencies 
(such as the median frequency, the zero crossing density, the so-called hi- 
5 over-low value, etc.). The recovery time constant depends mostly on the 
density of capillaries in the muscle and is rather insensitive to the exerted 
force. The fatigue time constant is strongly dependent on the force when it 
exceeds a certain critical level F c - The relation is: 

10 T F = n/(F-F c ) forF>F c (18a) 

and 

T F -> oo for F<F C (18b) 

15 

The combination of repeated work and recovery events causes the 
center frequency to decrease from the initial value to a final plateau value OF*,, 
at which there is a balance between the metabolite production during work and 
wash-out during recovery. The plateau value is: 

20 

CF«, = CF 0 (1 -k)T f /[(1 -k)T f + kTr] (19) 
Introducing the notations: 
25 ACF = CFo-CF« (20) 

and 

E = ACF/CFo (21) 

30 

Equation (19) can then be rearranged to read: 
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k = 1 / [1 + (T R / T F )(ACF / CFo. ) ] (22) 
With the notation: 

5 

Q = T R /T F (23) 
it is found that: 

10 Q = [(1 -k)/k] ACF/CFco (24) 

which is an experimentally measurable quantity. 
Force and pressure 



15 
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The force F can be determined for skeletal muscles working over joints 
without synergistic effects from other muscles. For the diaphragm muscle the 
force cannot be directly measured, rather the transdiaphragmatic pressure P di 
is obtained as a proportional measure. The following relation could be used: 

F = ^ E (25) 



where \i is a proportionality constant and E is the myoelectric signal strength, 
preferably based on the first spectral moment which is rather insensitive to 
25 metabolic changes caused by fatigue. The relation to the pressure is 
proportional but non-linear. This fact is taken into consideration by introducing 
the factor G(V) which is volume (V) dependent, i.e.: 

Pdi = FG(V) (26) 

30 

Thus, 



\ 
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fiG(V) = P dj /E (27) 

which also is an experimentally measurable quantity. 

Myoelectric signal strength and spectral changes 

Rearrangement of equation (18a) and insertion of equations (23) and 
(25) gives: 



a E - P - Q = 0 (28) 

where 

15 a = *iT r /ti (29) 

and 

P = FcTr/ti (30) 

20 

It can be observed that a is dependent (through the parameter jx) on 
the electrode geometry and placement in relation to the muscle, while the 
other parameters are rather constant for similar muscles. 

25 . Experiments under fatiguing conditions at any volume give 

corresponding values of E and Q (through the center frequency changes). A 
data fitting procedure (not regression) gives numerical values to a and p. With 
a and p known, an estimate of the muscle force F can be obtained^n relation 
to Its fatigue threshold value, i.e.: 



F/F c = Ea/p (31) 
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As long as F/F c is smaller than one, isometric fatigue does not develop. 
That means that the signal strength should be lower than the critical value: 



For periodic muscle work, higher forces and signal levels are tolerable. 

Spectral changes as indicators of tolerable concentration of metabolites 

The relative spectral change e, defined in equation (21), is an indirect 
measure of remaining concentration of metabolites in the muscle during 
periodic fatiguing contractions. It seems that the muscle very rapidly goes into 
an anaerobic metabolic state once the force is higher than F c and that virtually 
all contractions above this level causes changes of the center frequency. 
Therefore it is likely that a certain small value of e is tolerable as long as it is 
below a certain critical level, which we denote e c . With this critical value 
introduced into equation (22) and simultaneous use of equations (18a) and 
(23), it can be found that the condition for long term fatigue not occur to be: 



E < Ejsom = (3/cc 



(32) 



k < 1/{1+[(1 -ec)/e c ]T R (F-Fc)/r,} 



(33) 



This expression can be rearranged to give the force condition: 



F < F c +[(1-K)/K][e c /(1-e c )]Ti/T R 



(34) 



or, together with equation (30), 



F < F c {1 + [(1-K)/K][6 C /(1~e c )]/0} 



(35) 



Since the force in diaphragmatic contractions cannot be simply measured, 
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equations (33) to (35) are expressed as functions of the myoelectric signal 
strength and the transdiaphragmatic pressure. Use of equations (25) and (26) 
give for the myoelectric signal strength: 

k < 1/{1+[(1 -ec)/e c ](<xE-P)} (36) 

and 

E < {p + [(1 -k)/k][6 C /(1 -e c )]}/a (37) 
and for the transdiaphragmatic pressure: 

k < 1/{1 +[(1 -ec)/ec](aP di -0)} (38) 

and 

Pdi < |iG(V){P + [(1-k)/k][8c/(1-cc)]} (39) 
Procedures 

From equation (28) the coefficients a and p are determined with 
myoelectric data from fatigue tests (calibration). Fatigue test can be performed 
by either reducing the level of assist, or performing a short airway occlusion, 
while measuring the myoelectric activity during a few inspiratory attempts. To 
shorten and facilitate the fatigue test the subject could be encouraged to 
voluntarily increase his efforts. Such a test is routinely performed to determine 
the maximum inspiratory airway pressure. 

If myoelectric monitoring is used (giving signal strength and duty cycle), 
the critical signal strength is calculated using equation (37). 
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If transdiaphragmatic pressure monitoring is used (giving pressure, 
flow/volume, and duty cycle), the critical pressure level is calculated by using 
equation (39). This requires knowledge about the geometrical G(V) 
dependence, which can be obtained from a calibration of the experimentally 
5 measurable quantity Pdi/E as shown in equation (27). 

This is performed by either continuously measuring the P di through, for 
example, esophageal and gastric pressure sensors, and the volume through, 
for example, a flow meter. Geometrical dependence of inspiratory pressure 
10 can also be estimated by performing single or multiple breath airway 
occlusions at two lung volumes, e.g. end-inspiration and end expiration lung 
volumes, when the volume difference is measured with, for example, at least 
one flow meter, 

15 Estimates of the critical level £ C can be determined from the general 

experimental fact that fatigue does not occur below a duty cycle of 0.2 even at 
maximum muscle force and that the critical force level F c is approximately 0.2 
times the maximum force. Equation (33) then gives e c » p/(0+8/9), or, since 
both e c and 0 are small quantities: 

20 

ec « 0 (40) 

The gain of ventilatory assist should be adjusted at a level such that the 
myoelectric activity does not exceed that described in equation (37) (higher 
25 support suggest unnecessary muscle inactivation). However, myoelectric 
activity should not exceed that described in equation (32) (this level indicates 
the level for muscle fatigue during isometric contractions). 

Although the present invention has been described in relation to non- 
30 restrictive illustrative embodiments thereof, it should be kept in mind that these 
embodiments can be modified without departing from the spirit of the present 
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invention. In particular but not exclusively, the present invention pertains not 
only to CFdi and RMS but also to other related measures. Similarly, the 
present invention is concerned with any method of mechanical ventilation, 
including negative pressure ventilation. 

5 

REFERENCES 



[1] Aldrich TK, Sinderby C, McKenzie DK, Estenne M, and Gandevia SC. 
Electrophysiologic Techniques for the Assessment of Respiratory 
10 Muscle Function. In ATS/ERS Statement on respiratory muscle testing. 

.Am J Respir Crit Care Med 166: 610-623, 2002. 
[2] Bellemare F, and Grassino A. Effect of pressure and timing of 
contraction on human diaphragm fatigue. J Appl Physiol: Respirat 
Environ Exercise Physiol 53: 1 1 90-1 1 95, 1 982. 
15 [3] Bellemare F, and Grassino A. Evaluation of human diaphragm fatigue. J 
Appl Physiol: Respirat Environ Exercise Physiol 53: 1196-1206, 1982. 
[4] Beck J, Sinderby C, Lindstrom L, and Grassino A. Influence of bipolar 
electrode positioning on measurements of human crural diaphragm 
EMG. J Appl Physiol 81: 1434-1449, 1996. 
20 [5] Beck J, Sinderby C, Lindstrom L, and Grassino A. Diaphragm 
interference pattern EMG and compound muscle action potentials: 
effects of chest wall configuration. J. Appl. Physiol. 82: 520-530, 1997. 
[6] Beck J, Sinderby C, Lindstrom L, and Grassino A. Effects of. lung 
volume on diaphragm EMG signal strength during voluntary 
25 contractions. J Appl Physiol 85: 1 123-1 134, 1998. 

[7] Beck J, Sinderby C, Weinberg J, and Grassino A. Effects of muscle-to- 
electrode distance on the human diaphragm electromyogram. J Appl 
Physiol 79: 975-985, 1 995. 
[8] Brochard L, Harf A, Lorino H, and Lemaire F. Inspiratory pressure 
30 support prevents diaphragmatic fatigue during weaning from mechanical 

ventilation. Am Rev Respir Dis 139: 513-521, 1989, 



Copy provided by USPTO from the IFW Image Database on 01/31/2005 



[9] Broman, H. An investigation on the influence of a sustained contraction 

on the succession of action potentials from a single motor unit. 

Electromyogr Clin Neurophysiol 17:341-58, 1977. 
[10] Calzia E, Lindner KH, Witt S, Schirmer U, Lange H, Stenz R, and 
5 Georgieff M. Pressure-time product and work of breathing during 

biphasic continuous positive airway pressure and assisted spontaneous 

breathing. Am J Respir Crit Care Med 150: 904-910, 1994. 
[11] Clanton TL, Hartman E, and Julian MW. Preservation of sustainable 

inspiratory muscle pressure at increased end-expiratory lung volume. 
1 0 Am Rev Respir Dis 1 47: 385-391 , 1 993. 

[12] Clausen T, and Everts ME. K+ induced inhibition of contractile force in 

rat skeletal muscle, role of Na+ -K+ transport. Am J Physiol 261 (Cell 

Physiol. 30): C799-C807, 1991. 
[13] Cohen CA, Zagelbaum G, Gross D, Roussos C, and Macklem PT. 
15 Clinical manifestations of respiratory muscle fatigue. Am J Med 73: 308- 

316,1982, 

[14] Farkas GA and Roussos CH. Acute diaphragmatic shortening: In vitro 

mechanics and fatigue. Am Rev Respir Dis 130: 434-438, 1984. 
[15] Gandevia SC, and McKenzie DK. Human diaphragmatic EMG: changes 
20 with lung volume and posture during supramaximal phrenic nerve 

stimulation. J Appl Physiol 60: 1420-1428, 1986. 
[16] Gross D, Grassino A, Ross WRD, and Macklem PT. Electromyogram 

pattern of diaphragmatic fatigue. J Appl Physiol: Respirat Environ 

Exercise Physiol 46: 1-7, 1979. 
25 [17] Hodgkin AL. A note on conduction velocity. J Physiol (Lond) 125: 221- 

224, 1954. 

[18] Hodgkin AL and Huxley AF. A quantitative description of membrane 
current and its application to conduction and excitation in nerve. J 
Physiol (Lond) 117: 500-544, 1952. 
30 [19] Hussain S. Regulation of ventilatory muscle blood flow. J Appl Physiol 
81: 1455-1468, 1996. 



\ 

Copy provided by USPTO from the IFW Image Database on 01/31/2005 



28 



[20] Jubran A, Van de Graaff WB, and Tobin MJ. Variability of patient- 
ventilator interaction with pressure support ventilation in patients with 
chronic obstructive pulmonary disease. Am J Respir Crit Care Med 152: 
129-136, 1995. 

5 [21] Kadefors R, Kaiser E, and Petersen I. Dynamic spectrum analysis of 
myo-potential with special reference to muscle fatigue. Electromyog Clin 
Neurophysiol 8: 39-74, 1968. 
[22] Klawitter PF, and Clanton TL. Tension-time index, fatigue, and 
energetics in isolated rat diaphragm: a new experimental model. J Appl 
1 0 Physiol 96: 89-95, 2003. 

[23] Korner L, Parker P, Almstrom C, Herberts P, and Kadefors R. The 
relation between spectral changes of the myoelectric signal and the 
intramuscular pressure of the human skeletal muscle. Eur J Appl Physiol 
52: 202-206, 1984. 

15 [24] Laghi F, Cattapan SE, Jubran A, Parthasarathy S, Warshawsky P, Choi 
Y-S A, and Tobin MJ. Is weaning failure caused by low- frequency 
fatigue of the diaphragm? Am J Respir Crit Care Med 167: 120-127, 
2003. 

[25] Lindinger Ml, and Sjo gaard G. Potassium regulation during exercise 
20 and recovery. Sports Med 1 1 :382-401 , 1 991 . 

[26] Lindstrom, L. Fatigue changes in the myoelectric signal during periodic 

muscle work. Bull Eur Physiopathol Respir 15 Suppl: 107-114, 1979. 
[27] Lindstrom, L and Hellsing, G. Masseter muscle fatigue in man 
objectively quantified by analysis of myoelectric signals. Arch Oral Biol 
25 28:297-301,1983. 

[28] Lindstrom L, Kadefors R, and Petersen I. An electromyographic index 
for localized muscle fatigue. J Appl Physiol: Respirat Environ Exercise 
Physiol 43: 750-754, 1 977. 
[29] Lindstrom L and Magnusson R. Interpretation of myoelectric power 
30 density spectra: a model and its application. Proc IEEE 65: 653-662, 

1977. 



Copy provided by USPTO from the IFW Image Database on 01/31/2005 



29 



[30] Lindstrom L, and Petersen I. Power spectrum analysis of EMG signals 
and its applications. In: Progress in Clinical Neurophysiology. Computer- 
Aided Electromyography, edited by Desmedt JE. Basel: Karger, vol. 10, 
1983 p. 1-51. 

5 [31] Mortimer JT, Magnusson R, and Petersen I. Conduction velocity in 
ischemic muscle: effect on EMG frequency spectrum. Am J Physiol 219: 
1324-1329, 1970. 

[32] Ranieri VM, Giuliani R, Mascia L, Grasso S, Petruzzelli V, Puntillo N, 
Perchiazzi G, Fiore T, and Brienza A. Patient-ventilator interaction 
10 during acute hypercapnia: pressure-support vs. proportional-assist 

ventilation. J Appl Physiol 81:426-36, 1996. 
[33] Roussos C, Fixley M, Gross D, and Macklem PT. Fatigue of inspiratory 
muscles and their synergistic behavior. J Appl Physiol: Respirat Environ 
Exercise Physiol 46: 897-905, 1979. 
15 [34] Roussos CS and Macklem PT. Diaphragmatic fatigue in man. J Appl 
Physiol: Respirat Environ Exercise Physiol 43: 189-197, 1977. 
[35] Sasson CSH, Light RW, Lodio R, Siek GC, and Mahutte CK. Pressure- 
time product during continuous positive airway pressure, pressure 
support ventilation, and T-piece during weaning from mechanical 
20 ventilation Am Rev Respir Dis 143: 469-475, 1991. 

[36] Sinderby C, Beck JC, Lindstrom L, and Grassino A. Enhancement of 
signal quality in esophageal recordings of diaphragm EMG. J Appl 
Physiol 82: 1370-1377, 1997. 
[37] Sinderby C, Beck J, Weinberg J, Spahija J, and Grassino A. Voluntary 
25 activation of the human diaphragm in health and disease. J Appl Physiol 

85:2146-2158, 1998. 
[38] Sinderby CA, Comtois AS, Thomson RG, and Grassino AE. Influence of 
the bipolar electrode transfer function on the electromyogram power 
spectrum. Muscle & Nerve 19: 290-301, 1996. 
30 [39] Sinderby C, Lindstrom L, Comtois N, and Grassino AE. Effects of 
diaphragm shortening on the mean action potential conduction velocity 



Copy provided by USPTO from the IFW Image Database on 01/31/2005 



30 



in canines. J Physiol 490: 207-214, 1996. 
[40] Sinderby C, Lindstrom L, and Grassino A. Automatic assessment of 

electromyogram quality. J Appl Physiol 79: 1803-1815, 1995. 
[41] Sinderby C, Navalesi P, Beck J, Skrobik Y, Comtois N, Friberg S, 
5 Gottfried SB, and Lindstrom L. Neural control of mechanical ventilation. 

Nat Med 5: 1433-1436, 1999. 
[42] Sinderby C, Spahija J, and Beck J. Changes in respiratory effort 

sensation over time are linked to the frequency content of diaphragm 

electrical activity. Am J Respir Crit Care Med 163: 905-910, 2001 . 
10 [43] Tobin MJ, Perez W, Guenther SM, Semmes BJ, Mador MJ, Allen SJ, 

Lodato RF, Dantzker DR. The pattern of breathing during successful 

and unsuccessful trials of weaning from mechanical ventilation. Am Rev 

Respir Dis 1 34: 1 1 1 1 -1 1 1 8, 1 986. 
[44] Tzelepis G, McCool FD, Leith DE, and Hoppin FG Jr. Increased lung 
15 volume limits endurance of inspirator y muscles. J Appl Physiol 64: 

1796-1802, 1988. 



Copy provided by USPTO from the IFW Image Database on 01/31/2005 



Table 1. Breathing pattern and targeted Pdi values during the three maneuvers performed. 



Volume maneuver 
Subject Ti/Ttot VT%IC Pdi 



Expulsive maneuver Expulsive maneuver 
Low pressure High pressure 



Ti/Ttot Pdi 



Ti/Ttot Pdi 



0.66 



89.2 



47.8 



0.66 



45.1 



0.67 



77.3 



0.66 



65.6 



25.3 



0.67 



26.6 



0.66 



93.3 



0.67 



68.3 



10.2 



0.65 



11.6 



0.66 



39.3 



0.67 



88.8 



38.6 



0.66 



39.8 



0.64 



77.3 



Mean 
(±SD) 



0.65 

0.66 
(0.01) 



75.5 

77.5 
(11.1) 



34.1 

31.2 
(14.2) 



0.67 

0.66 
(0.01) 



35.0 

31.6 
(13.1) 



0.67 

0.66 
(0.01) 



50.3 

67.5 
(22.0) 



Values are means for each subject of all the maneuvers performed. Ti/Ttot, duty cycle; Pdi, 
transdiaphragmatic pressure, Vr, tidal volume; IC, inspiratory capacity. All subjects were able to 
maintain the imposed duty cycle. 
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Table 2. Individual CFdi values observed at end of each maneuver. 



Subject 


Volume maneuver 


Expulsive maneuver 
Low pressure 


Expulsive maneuver 
High pressure 


1 


68.3 ± 8.7 


81.3 ±7.9 


61.6 ±5.6 


2 


67.4 ±4.5 


86.7 ±8.0 


70.5 ± 13.2 


3 


82.0 ± 12.3 


94.3 ±7.8 


71.5 ±5.0 


4 


72.8 ±5.8 


89.4 ±6.1 


73.4 ±3.1 


5 


80.4 ±9.3 


104.1 ±4.4 


83.6 ± 5.4 



Mean±SD 74.2±6.8 91.2±8.6 72.1+7.9 



Values are means for each subject for each of the maneuvers performed. 
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Table 3. Test-retest of the volume maneuver 



CFdio CFdioo EAdi Pdi 

(Hz) (Hz) (% max) (cm H 2 0) 



Subject 


Voll 


VoI2 


VoU 


VoE 


Voll 


Vol2 


Voll 


Vol2 


1 


90.4 


93.1 


68.3 


66.2 


46.5 


38.1 


47.8 


44.6 


2 


102.5 


108.2 


67.4 


65.7 


42.4 


44.4 


38.6 


38.3 


3 


94.7 


100.3 


82.0 


84.7 


49.8 


55.0 


34.1 


30.7 


4 


99.7 


99.5 


72.6 


72.1 


51.3 


50.6 


36.5 


35.8 


Mean 


96.8 


100.3 


72.8 


71.9 


66.7 


65.0 


25.3 


29.5 


(±SD) 


(5.4) 


(6.2) 


(6-7) 


(8.8) 


(10.7) 


(11.9) 


(9.4) 


(7.0) 


ICC 


0.94 


0.98 




0.93 




0.95 





EAdi, diaphragm electrical activity calculated as root-mean-square; CFdi 0 , baseline center 
frequency determined during resting breathing; CFdioo, plateau value of the center frequency at 
the end of the volume maneuver; Pdi, transdiaphragmatic pressure, Voll, first volume maneuver 
performed, Vo 12, second volume maneuver performed, ICC, interclass correlation coefficient. 
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Table 4. Experimental results and calculated values from Appendix L 



Force ratios Geometric factors 





Pdi EAdiCFo CF„ 
(cm H 2 0) (a.u.) (Hz) (Hz) 


Q 


«t> ■ <t>ni 


Yi Yu 
(cmH 2 0)(cm H z O) 


Ym 


Volume Maneuver 


31.2 


60.9 100 74.2 


0.175 


6.12 5.20 5.66 


122 


144 


133 


Expulsive maneuver 
(Low pressure) 


31.6 


24.9 100 91.2 


0.050 


2.50 2.19 2.35 


303 


347 


325 


Expulsive maneuver 
(High pressure) 


67.5 


44.3 100 72.1 


0.196 


4.45 5.69 5.07 


364 


285 


324 



Pdi, transdiaphragmatic pressure; a.u., arbitrary units; EAdi, electrical activity of the diaphragm; 
CFdi o, baseline diaphragm center frequency determined during resting breathing; CFdi ««, 



plateau value of the diaphragm center frequency at the end of the maneuver; Q, ratio of the time 
constants of CFdi recovery and decline, see equation (5); i, see equation (9); § n , see equation 
(10); ()> mi meanof<j>i and(|>n; Yh see equation (14); y ii> see equation (15);Ym, meanofyi andyn. 
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Volume maneuver (Low Pdi) 
Expulsive maneuver (Low Pdi) 
Expulsive maneuver (High Pdi) 
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* f*<0.05 compared 
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